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Abstract 

Wolbachia is the most prevalent symbiont described in arthropods to date. Wol- 
bachia can manipulate host reproduction, provide nutrition to insect hosts and 
protect insect hosts from pathogenic viruses. So far, 13 supergroups of Wolba- 
chia have been identified. The whitefly Bemisia tabaci is a complex containing 
more than 28 morphologically indistinguishable cryptic species. Some cryptic 
species of this complex are invasive. In this study, we report a comprehensive 
survey of Wolbachia in B. tabaci and its relative B. afer from 1658 insects repre- 
senting 54 populations across 13 provinces of China and one state of Australia. 
Based on the results of PGR or sequencing of the 16S rRNA gene, the overall 
rates of Wolbachia infection were 79.6% and 0.96% in the indigenous and inva- 
sive Bemisia whiteflies, respectively. We detected a new Wolbachia supergroup 
by sequencing five molecular marker genes including 16S rRNA, groEL, gltA, 
hcpA, and fbpA genes. Data showed that many protein-coding genes have limi- 
tations in detecting and classifying newly identified Wolbachia supergroups and 
thus raise a challenge to the known Wolbachia MLST standard analysis system. 
Besides, the other Wolbachia strains detected from whiteflies were clustered into 
supergroup B. Phylogenetic trees of whitefly mitochondrial cytochrome oxidase 
subunit I and Wolbachia multiple sequencing typing genes were not congruent. 
In addition, Wolbachia was also detected outside the special bacteriocytes in 
two cryptic species by fluorescence in situ hybridization, indicating the horizon- 
tal transmission of Wolbachia. Our results indicate that members of Wolbachia 
are far from well explored. 



Introduction 

Wolbachia are rickettsial endosymbiotic bacteria in the 
class Alphaproteobacteria. Wolbachia bacteria are consid- 
ered the most widespread endosymbionts in animals as 
they are found in all major classes of arthropods and 
some nematodes (Jeyaprakash and Hoy 2000; Werren and 
Windsor 2000; Duron et al. 2008; RusseU et al 2012). A 
meta-analysis suggests that the proportion of Wolbachia 
infection in insect species in the terrestrial world is about 
40% (Zug and Hammerstein 2012). 

In some host species, the successful maintenance and 
spread of Wolbachia is mainly achieved by the induction 
of cytoplasmic incompatibility to produce more female 
offspring, thus enhancing its maternal transmission 
(Stouthamer et al. 1999). In addition, manipulation of 
reproduction by Wolbachia includes feminizing genetic 
males, causing parthenogenesis, and killing male progenies 



(Stouthamer et al. 1999; Werren et al. 2008). Recent stud- 
ies found that Wolbachia benefits insect hosts by providing 
essential nutrition (Hosokawa et al. 2010), enhancing host 
stem cell's proliferation (Fast et al. 2011), and protecting 
insect from pathogenic RNA viruses (Hedges et al. 2008). 

The genus Wolbachia is highly divergent and has so far 
been divided into 13 supergroups (A-N, except for G 
which is a combination of A and B) (Lo et al. 2002, 2007; 
Baldo and Werren 2007; Haegeman et al. 2009; Ros et al. 
2009; Augustinos et al. 2011). Wolbachia supergroups are 
characterized mainly with molecular markers such as rrs 
(I6S rRNA), ftsZ (cell division protein), gltA (Gitrate syn- 
thase), groEL (Ghaperonin GroEL) and wsp (Wolbachia 
surface protein) genes (O'Neill et al. 1992; Zhou et al. 
1998; Werren and Windsor 2000; Gasiraghi et al. 2005). 
Wolbachia genotyping is inferred mainly from multi locus 
sequence typing (MLST) genes (gatB, coxA, hcpA, fbpA, 
and ftsZ genes) and amino acid sequences of the four 
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hypervariable regions (HVRs) of WSP protein (Baldo 
et al. 2005, 2006). 

Bemisia tabaci (Hemiptera: Aleyrodidae) is a complex 
containing more than 28 morphologically indistinguish- 
able cryptic species (De Barro et al. 2011; Hu et al. 2011). 
Through millions of years of evolution, the various cryptic 
species of this complex show a clear geographic pattern of 
distribution around the globe (Boykin et al. 2007, 2013; 
De Barro et al. 2011). However, with the development of 
modern transport, whiteflies have been transferred fre- 
quently among different continents (Naranjo et al. 2010). 
During the last twenty years, two cryptic species of the 
B. tabaci complex, Middle East-Asia Minor 1 (formerly 
known as the B "biotype," hereafter MEAMl) and Medi- 
terranean (formerly known as the Q "biotype," hereafter 
MED) have invaded many regions of the world (Dalton 
2006; Hu et al. 2011). They have caused serious damages 
to local agriculture through direct plant sap sucking and 
transmission of plant pathogenic viruses (Oliveira et al. 
2001). What is more, the rapid invasion of MEAMl and 
MED has caused the replacement of many indigenous 
cryptic species of the B. tabaci complex (Liu et al. 2007; 
Hu et al. 2011; Muniz et al. 2011; Rao et al. 2011). These 
events provide us a unique opportunity for studying the 
evolution and transmission of Wolbachia among different 
B. tabaci cryptic species, which were geographically iso- 
lated in history but have become sympatric recently. 

Previous studies have investigated the diversity of Wol- 
bachia in the B. tabaci species complex (Nirgianaki et al. 
2003; Chiel et al. 2007; Gueguen et al. 2010; Chu et al. 
2011; Pan et al. 2012; Singh et al. 2012; Bing et al. 2013a). 
However, most of these reports focused on the two invasive 
cryptic species MEAMl and MED and only used one to 
three marker genes in the investigation, and the distribu- 
tion of Wolbachia in most indigenous B. tabaci cryptic spe- 
cies remains largely unknown. In this study, we examined 
the distribution of Wolbachia in B. afer and 10 cryptic spe- 
cies of the B. tabaci species complex collected from 13 
provinces of China and one state of Australia. We report: 

(1) the prevalence of Wolbachia in B. afer and B. tabaci; 

(2) the discovery of a probably new Wolbachia (supergroup 
O) in whiteflies by sequencing of rrs gene and four protein- 
coding genes {fbpA, hcpA, gltA, and groEL); (3) the diver- 
sity and phylogenetic status of Wolbachia strains within 
these whiteflies; and (4) evidence for horizontal transfer of 
Wolbachia among B. tabaci cryptic species. 

Materials and Methods 

Whitef ly collection and DNA extraction 

Bemisia specimens were collected from 13 provinces of 
China and one state of Australia. Details for collection 



(geographical locations, host plants, and dates) of those 
populations are summarized in Fig. 1 and Table Al. 
Whiteflies collected from the same locality and host plant 
were considered as one population. Whitefly samples were 
initially immersed in 95% ethanol after collection and 
subsequently kept at — 20°C until DNA extraction. Total 
whitefly DNA was extracted from individual adult speci- 
mens according to the method of DeBarro and Driver 
(1997). The quality of the DNA samples was confirmed 
by PCR amplification of a 0.8 kb fragment of whitefly 
mitochondrial cytochrome oxidase I (mtCOl) gene using 
the primers Cl-J-2195 and L2-N-3014 (Table A2). Cryp- 
tic species of B. tabaci were first identified based on the 
polymerase chain reaction restriction fragment length 
polymorphism (PCR-RFLP) method described by Qin 
et al. (2013), and the sex of whiteflies was identified 
through genital morphology. A total of 1658 whitefly 
DNA samples were positive for PCR amplification using 
the mtCOI primers, indicating satisfactory quality of the 
DNA templates. 

Diagnostic screening of Wolbachia 

The presence of Wolbachia was screened based on the 
amplification of a 0.6 kb fragment with the Wolbachia rrs 
primers (Table A2). Standard PCR analyses were per- 
formed using 2xEasyTaq PCR SuperMix (TransGen, 
Beijing, China) in a PTC-200 Thermocycler (Bio-Rad, 
Hercules, CA). PCR procedures were an initial step of 
94°C for 3 min, followed by 32 cycles of 94°C for 30 s, 
55°C for 30 s, and 72°C for 45 s and a final step of 72°C 
for 10 min. Amplified DNA products were electrophoresed 
on agarose gels and stained with GelRed (Biotium, San 
Francisco, CA). To verify PCR results, amplified bands 
(especially uncertain ones) were purified by AxyPrep DNA 
gel extraction kit (Axygen, Silicon Valley, CA) and cloned 
into the pGEM-T vector (Promega, Madison, WI). Plas- 
mids containing the DNA inserts of expected sizes were 
confirmed by PCR and sequenced in an ABI 3730 DNA 
analyzer (Applied Biosystems, Foster City, CA). Sequenc- 
ing results were then checked by Blast in NCBI nr database 
(http://blast.ncbi.nlm.nih.gov/Blast.cgi). Only those indi- 
viduals, which were blasted to expected products of the 
specific primers, were considered to be infected. All PCRs 
included a negative control (sterile water instead of DNA) 
and a positive control (DNA of China 1 whitefly). For 
Wolbachia-positive populations, Bemisia species were fur- 
ther identified by phylogenetic analysis of the mtCOI gene 
(Fig. A3). Wolbachia infection rates between whitefly sexes 
were statistically tested using the Fisher's exact test. Statis- 
tical significance of the infection rates among different 
B. afer and B. tabaci cryptic species was calculated using 
the test and corrected by the Bonferroni procedure. AU 
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Figure 1. Localities of sampling and infection 
frequencies of Wolbachia in 53 Chinese 
populations of Bemisia. About 30 individuals 
from each population were subjected to 
diagnostic PCR analysis. Whiteflies collected 
from the same host plant in the same locality 
were considered as one population. The Arabic 
numerals correspond to populations numbered 
in Table Al. Figures in parentheses indicate 
the numbers of individuals sampled from each 
of the populations. Different colors represent 
6. afer and different cryptic species of the 
6. tabaci complex. The "#" signs indicate the 
laboratory lines that had been maintained on 
cotton since collection, and the "*" signs 
indicate the 5 populations that are positive for 
Wolbachia supergroup O (referred to in 
Table 3 and Fig. 2). 



Statistical analyses were performed using the Data Process- 
ing System (DPS) software (Tang and Zhang 2013). 

Sequencing and typing of Wolbachia 

As the phylogenetic analysis of the 0.6 kb rrs sequences 
indicated a possible new supergroup of Wolbachia, we 
amplified the rrs gene from whitefly populations and 
introduced an endonuclease Vspl (AT/TAAT) (Fisher 
Scientific, Pittsburgh, PA) to digest the target bands, to 
investigate the infection prevalence of the new supergroup 
of Wolbachia. Several whitefly individuals were then ran- 
domly selected for sequencing confirmation of PCR-RFLP 
results. The groEL, wsp, and MLST {gatB, coxA, hcpA, 
fbpA, and ftsZ) genes from every combination of host 
species and rrs genotype were amplified by TransTaq-T 
DNA Polymerase (TransGen), cloned into pEASY-Tl vec- 
tors (TransGen) and sequenced on ABI 3730 DNA ana- 
lyzer (Applied Biosystems). 

In addition, to confirm the finding of the new Wolba- 
chia supergroup, nearly a complete rrs sequence 
(1417 bp) and a partial gltA sequence (659 bp) were 
amplified from one Wolbachia new supergroup singular- 
infected population. PCR amplifications, DNA cloning 
and sequencing procedures were accomplished as 
described previously. The cycling procedures were the 
same as described earlier with changes on the annealing 
temperature for different primers. Primer sequences and 
annealing temperatures of rrs, gltA, groEL, wsp, and the 
MLST genes were listed in Table A2. 



Sequences of MLST and wsp genes were manually 
trimmed in line with the template provided in Wolbachia 
MLST website and compared with sequences in the Wol- 
bachia MLST database (http://pubmlst.org/wolbachia/). 
Novel sequences were submitted to the database curators 
as new alleles. Each unique combination of five MLST 
sequences was designated a strain type (ST) number in 
the Wolbachia MLST database (Baldo et al. 2006). Previ- 
ously published sequences from other whiteflies were 
added to the data set to increase the power of statistical 
comparisons. All newly obtained allele numbers and ST 
numbers in this study are summarized in Table 1. 

iVIolecular phylogenetic analysis 

Phylogenetic analyses were constructed using (1) Wolba- 
chia sequences of rrs, gltA, groEL, MLST, and wsp 
genes of different supergroups from various hosts; and 
(2) whitefly mtCOI sequences. All sequences used in this 
study were edited and aligned manually using Clustal W 
(ver. 1.6) (Thompson et al. 1994) in MEGA (ver. 5.10) 
(Tamura et al. 2011). The Gblocks program (ver. 0.91b) 
(Castresana 2000) was used to remove poorly aligned 
positions and to obtain nonambiguous sequence align- 
ments. The best-fit evolutionary model for the sequence 
data was determined using hierarchical likelihood ratio 
tests and Akaike information criterion with the program 
jModelTest (ver. 0.1.1) (Posada 2008). Phylogenetic trees 
were constructed with the Bayesian inference using 
MrBayes (ver 3.1.2) (Ronquist and Huelsenbeck 2003). 
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For these gene data, 5 million generations were run; 
50,000 trees were obtained, and the first 25% trees were 
discarded as burn-in. The resulting phylogenetic trees 
were visualized in TreeView (ver. 1.6.6) (Page 1996). The 
comparison between the phylogeny of whitefly mtCOI 
and Wolbachia concatenate MLST data sets was con- 
structed with Dendroscope (ver. 3.2.8) (Huson and 
Scornavacca 2012). The new Wolbachia rrs, gltA, groEL, 
wsp, and MLST gene sequences and whitefly mtCOI gene 
sequences have been deposited in the GenBank database 
(Table 1, Fig. A3 and Table A4). 

The pairwise genetic divergence of different Wolbachia 
supergroups was calculated using the Kimura 2-parameter 
method (Kimura 1980) in MEGA (ver. 5.10) (Tamura 
et al. 2011). Because recombination of sequences has 
potentially disruptive influences on phylogenetic-based 
molecular evolution analyses (Martin et al. 2011), align- 
ments of individual and concatenated genes were checked 
for significant levels of recombination using the Phi test 
(Bruen et al. 2006) in SplitsTree4 under default condi- 
tions (Huson and Bryant 2006). When recombination 
was tested to be significant, a phylogenetic network 
framework was constructed based on uncorrected P dis- 
tances using the Neighbor-net method (Bryant and Moul- 
ton 2004) implemented in SplitsTree4 (ver. 4.13.1) 
(Huson and Bryant 2006). 

FISH 

Localization of Portiera and Wolbachia was studied in 
nymphs and adults of B. tabaci Asia II 1 (Pop. 10) and 
Asia II 9 (Pop. 43) using fluorescence labeled probes spe- 
cifically targeting the rrs genes of these bacteria. We fol- 
lowed the previous protocols for the FISH experiments 
(Bing et al. 2013b). Briefly, specimens were collected 
directly into Carnoy's fixative and fixed overnight. After 
fixation, the samples were hybridized overnight in hybrid- 
ization buffer (20 mM Tris-HCl, pH 8.0, 0.9 M NaCl, 
0.01% sodium dodecyl sulfate, 30% deionized formam- 
ide) containing 10 pmol of fluorescent probes. The probe, 
BTP1-Cy3 (5'-Cy3- TGTCAGTGTCAGCCCAGAAG-3'), 
was used to target rrs gene of Portiera (Gottlieb et al. 
2006). A new probe, Wolb-1-488 (5'- Alexa Fluor 488- 
TAATATAGGCTCATCTAATAGCAA -3'), was designed 
to target rrs gene of Wolbachia. The specificity of the 
detection was first checked by "probe match" in RDP 10 
(update to May 14, 2013) (Cole et al. 2009) and BLAST 
in nr database of NCBI (http://blast.ncbi.nlm.nih.gov/ 
Blast.cgi) and then confirmed using the following con- 
trols: a no probe control and Wolbachia-free whiteflies 
(samples of the B. tabaci MED species and MEAMl spe- 
cies). Stained samples were wholly mounted and viewed 
under a Zeiss LSM780 confocal microscope. 



Results 

Prevalence of Wolbachia in Bemisia species 

Bemisia tabaci has a wide distribution in China. In this 
study, samples of Bemisia were obtained from 24 localities 
of 13 provinces of this country (Fig. 1). In all, these sam- 
ples represent one population of B. afer and 52 popula- 
tions of 9 cryptic species of the B. tabaci complex from 
China (Table Al). In addition, one sample of an indige- 
nous species of the B. tabaci complex was obtained from 
Australia (Table Al). 

Of the 1658 individuals examined, rrs PCR assays indi- 
cated that the infection rates of Wolbachia varied among 
species, and even among populations of a given species, 
ranging from 0% to 100% (Fig. 1; Table 2, Chi-square 
test, P < 0.0001), but did not differ between sexes in each 
of the populations that were tested statistically 
(Table Al). The incidence of Wolbachia infection in 
indigenous whiteflies (79.61%, n = 618) was significantly 
higher than that in invasive whiteflies (MEAMl and 
MED, 0.96%, n = 1040; Fisher's exact 2-tailed test, 
P < 0.0001). The infection rate of Wolbachia also varied 
among different indigenous species of the B. tabaci com- 
plex. For instance, 88.4% of China 1 whiteflies were posi- 
tive for Wolbachia while so were only 2.1% of Asia II 3 
whiteflies (Table 2). The rate of Wolbachia infection in 
B. afer was 77.5% (Table 2). 

Diversity of Wolbachia infections 

For those Wolbachia-^ositvfe populations, 2-3 individuals 
were further analyzed by sequencing Wolbachia rrs gene 
(592 bp) and performing Bayesian phylogenetic analysis. 
Most of whitefly Wolbachia rrs sequences were clustered 



Table 2. Rates of Wolbachia infection in Bemisia afer and cryptic spe- 
cies of tine 6. tabaci complex. 



WInitefly 




No. of 


Infection 


species 


n 


localities 


rate (%)' 


6. afer 


40 


1 


77.5 a 


Asia 1 


123 


2 


99.2 b 


Asia II 1 


151 


4 


96.7 be 


Asia II 3 


48 


2 


2.1 de 


Asia II 6 


60 


2 


23.3 d 


Asia II 7 


38 


1 


97.4 ab 


Asia II 9 


35 


1 


88.6 ab 


China 1 


112 


3 


88.4 ac 


Australia 


11 


1 


100.0 ab 


MEAMl 


334 


12 


0.0 e 


MED 


706 


26 


1.4 e 



'Figures followed by different letters differ significantly (adjusted 
P < 0.05, using Bonferroni corrections). 
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Figure 2. Phylogeny of the Wolbachia 
identified from Bemisia afer and cryptic species 
of tine 6. tabaci complex based on bacterial rrs 
gene sequences (592 sites). Wolbachia strains 
are characterized by the names of their host 
species. The tree was constructed using a 
TPMIuf + G substitution model for Bayesian 
analysis. Bayesian posterior probabilities are 
shown on the branches. Sequences obtained in 
this study are shown in bold. The bar indicates 
a branch length of 0.1 substitutions/site. The 
sequence names and GenBank accession 
numbers are listed in Tables A4 and A6. 



Nasonia vitripennis 
(M84686) 
Bryobia praetiosa 
M Bemisia tabaci 



wBt_29-1 38 




into the supergroup B (Fig. 2). However, several Wolba- 
chia rrs sequences obtained from four Asia II 1 popula- 
tions (wBt_10, wBt_28-2, wBt_29-2, wBt_iO-2) and one 
MED population {wBt_2) formed a strict and robust 
monophyletic clade (Fig. 2). 

Identification of Wolbachia supergroup O 

Further Bayesian phylogenetic analysis on a nearly com- 
plete rrs sequence (1317 bp) of the strange Wolbachia 
(wBf_10) revealed that this Wolbachia differed widely from 
other known Wolbachia (Fig. A5). The divergence of rrs 
between wBt_W and the supergroup M is 2.52%, which is 
the smallest of that between wBt_W and all previously 
reported Wolbachia supergroups (A to N) (Table A6). In 
view of these apparent differences, we proposed to name 
these Wolbachia as supergroup O temporarily. 

Our results thus showed the presence of two genetically 
distant Wolbachia supergroups in whiteflies. To clarify the 
detailed infection prevalence of Wolbachia O in Bemisia 
whiteflies, PCR-RFLP was introduced to digest all positive 



rrs PGR products by the restriction enzymes Vspl. No 
Vspl restriction site was found in rrs sequences of Wolba- 
chia supergroup O, whereas rrs amplicons from Wolbachi- 
a supergroup B could be digested into multiple bands by 
Vspl (Fig. A7). Whiteflies in populations 2 (MED species) 
and 10 (Asia II 1 species) are infected singly by Wolbachia 
O, whereas populations 28, 29, and 30 (all three are Asia 
II 1 species) are infected by both Wolbachia O and Wol- 
bachia B (Table 3). 

At least one of the eight tested protein-coding genes 
{gltA, groEL, MLST [gatB, coxA, hcpA, fbpA, ftsZ), and 
wsp genes) was successfully amplified and sequenced for 
all Wolbachia-'mitcttA populations in this study. Both 
neighbor-net analysis of fbpA, gltA, hcpA gene and Bayes- 
ian interference of groEL gene supported the existence of 
the Wolbachia supergroup O (w5f_10) (Figs. 3 and 4, and 
Figs. AlO and A12). However, it should be noted that the 
results of phylogenetic analyses with different genes were 
not always consistent. In particular, analysis of fbpA, 
groEL, and hcpA gene clustered some strange Wolbachia, 
such as wBt_29-2 and wBt_30-2, which were identified as 
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Table 3. Infection frequencies of tine Wolbachia O in five populations of the Bemisia tabaci complex. 



Pop. no. 


Cryptic 
species 


n' 


% without 
Wolbachia 
infection 


Single infection 
0 


(%) 

B 


Double 
infection (%) 


2 


MED 


29 


93.1 


6.9 






10 


Asia II 1 


44 




100 






28 


Asia II 1 


43 


6.9 


14.0 


62.8 


16.3 


29 


Asia II 1 


35 


2.9 


5.7 


65.7 


25.7 


30 


Asia II 1 


29 


3.4 


24.1 


3.4 


69.0 



Number of whitefly individuals collected from the five populations shown with asterisks in Fig. 1 and Table A1. 
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Figure 3. Phylogenetic position of the 
Wolbachia identified from the putative species 
Asia II 1 of the Bemisia tabaci complex based 
on bacterial gItA gene sequences (636 sites) 
using the Neighbor-net method. Each edge (or 
a set of parallel edges) corresponds to a split in 
the data set and has a length equal to the 
weight of the split. The sequence obtained in 
this study is shown in bold. The bar indicates a 
branch length of 0.1 substitutions/site. The 
names and sequence GenBank accession 
numbers are listed in Table A4. 



ivBt 10 



\ 



Q Onchocerca gibsoni 
Dirofilaria repens 

o.as; 

J Dipetalonema gracile 



Nasonia vitripennis (AY714796) 
Bryobia praetiosa 
Bryobia sarothamni 
wBt_43 wBt_35 wBt_28-^ 
l*Bf_49-1 wBt_3S 
wBt_S4 wSf_39 
wBt 41 



B 
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l*St_30-2 




Bryobia sp. ^ 



Nasonia wfripenn/s (AY714812) A 
Zootermopsis nevadensis |-| 



Litomosoides 
sigmodontis 



Figure 4. Phylogenetic position of the 
Wolbachia identified from Bemisia afer and 
6. tabaci putative species based on bacterial 
graft gene sequences (491 sites). Wolbachia 
strains are characterized by the names of their 
host species. The tree was constructed using a 
GTR + G substitution model for Bayesian 
analysis. Bayesian posterior probabilities are 
shown on the branches. The sequence 
obtained in this study is shown in bold. The 
bar indicates a branch length of 0.1 
substitutions/site. The names and sequence 
GenBank accession numbers are listed in 
Table 1 and Table A4. 



O by rrs gene, into supergroup B (Fig. 4, and Figs. AlO 
and A12). 

Sixteen STs were identified in whiteflies from this 
study, and all of them are new to the MLST database 
(Table 1). Though efforts were made, some PCRs failed 
when amplifying the MLST and wsp genes from super- 
group O- infected whiteflies (Table 1). As a result. 



sequences from supergroup O-infected whiteflies were 
excluded from phylogenetic analysis of the concatenated 
MLST sequences. Respective Bayesian interference of sep- 
arate gatB, coxA, ftsZ, and wsp genes showed that all 
Wolbachia detected in whiteflies belonged to supergroup 
B (Fig. 5, and Figs. A8, A9 and All). Neighbor-net analy- 
sis clustered the majority of Wolbachia into supergroup B 
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1.00 



Figure 5. Phylogenetic position of tine 
Wolbachia identified from Bemisia afer and 
6. tabaci putative species based on bacterial 
W5p gene sequences (512 sites). Wolbachia 
strains are cinaracterized by tine names of their 
Inost species. The two Drosophiia wsp 
sequences are the outgroups. The tree was 
constructed using a TIMS + G substitution 
model for Bayesian analysis. Bayesian posterior 
probabilities are shown on the branches. The 
sequence obtained in this study is shown in 
bold. The bar indicates a branch length of 0.1 
substitutions/site. The names and sequence 
GenBanl< accession numbers are listed in 
parentheses and Table 1 . 



1.00 



0.1 
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Laodelphax striatellus (AF020080) 



wBaJ\ 
wBt 26 



1 .00 J wBt_27 
0.571. wBt 35 



except for wBf_10 (Figs. AlO and A12). The hcpA genes 
from wBt_10 and fbpA from wBf_10 and wBt_28-2 
formed a separate branch that differs distinctly from all 
known reference sequences. 

Co-divergence between the divergence of 
Wolbachia supergroup B and whitef ly 
species 

The codivergence of Bemisia and Wolbachia supergroup B 
was assessed by studying the sequences of partial mtCOI 
gene and Wolbachia MLST genes. For those Wolbachia 
identified from B. tabaci, very poor congruence was 
found between the phylogenies of tntCOI and 
concatenated MLST genes (Fig. 6). The topology of MLST 
tree differs obviously from that of mtCOI. Whiteflies 
belonging to the same cryptic species harbored distant 
Wolbachia strains. For instance, two populations {wBt_35 



and wBt_3i8) of Wolbachia identified from Asia II 6 are 
clustered in different phylogenetic groups. 

Localization of Wolbachia in Bemisia tabaci 

The FISH of bacteria revealed that Portiera was seen exclu- 
sively in the bacteriocytes of whiteflies. In the tested 
nymphs, Wolbachia was strictly located in the bacteriocytes 
among the abundant Portiera (Fig. 7). Nevertheless, in the 
adults, Wolbachia was detected both outside and inside the 
bacteriocytes (Fig. 7). Signals of Wolbachia shown at the 
anterior pole of the oocytes of female adults indicate its 
vertical transmission (arrows marked in Fig. 7 D & H). 

Discussion 

Wolbachia is widely distributed among invertebrates and 
is considered as the most prevalent symbiont identified so 
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Figure 6. Comparisons of Bemisia and 
Wolbacliia phylogenies. A, the whitefly 
phylogeny constructed based on Bayesian 
analysis of mtCOl sequences (817 bp) as 
sinown in Fig. A3 using TIIV13 + I + G model. B, 
the Wolbachia phylogeny constructed based 
on Bayesian analysis of concatenated 
sequences of MLST genes (2079 bp) as shown 
in Table 1 using GTR + I + G model. Bayesian 
posterior probabilities are shown on the 
branches. Dashed lines connect hosts to their 
respective Wolbachia strains. The scale bar is in 
units of substitutions/site. 



far. Though several research groups have investigated the 
prevalence of Wolbachia in some cryptic species of the 
B. tabaci complex, our study represents the first compre- 
hensive analysis of Wolbachia infection among both inva- 
sive and indigenous cryptic species of the B. tabaci 
complex in Asia. In addition, compared with previous 
investigations, we used five more molecular markers in 
our analyses. 

Prevalence of Wolbachia varies between 
invasive and indigenous whiteflies 

In this study, Wolbachia infection rates in five (Asia I, 
Asia II 1, Asia II 7, Asia II 9, and China 1) of the seven 
Chinese indigenous species reached over 70%. In contrast, 
Wolbachia infection rate in the MED populations from 
China was only 1.4% (10/706), and no infection (0/334) 
was detected in all MEAMl populations from this coun- 
try. The low rates of Wolbachia infection in MEAMl and 
MED agree with those observed in most previous studies. 
For example, in populations of MEAMl and MED from 
Europe and Western Africa, infection rates of Wolbachia 
varied from 0-8.3% and 0-33% (Nirgianaki et al. 2003; 
Chiel et al. 2007; Gueguen et al. 2010; Skaljac et al. 2010; 
Chu et al. 2011; Thierry et al. 2011; GnankinE et al. 
2013). And in populations of MEAMl and MED from 
China, the rates of Wolbachia infection were 0.2% (1/456) 
and 0% (0/1149), respectively (Pan et al. 2012). As a 
whole, our data indicate a high variability of prevalence 
of Wolbachia between cryptic species of the B. tabaci 
complex. In our sampling, we obtained adequate numbers 
of whitefly individuals for five (Asia II 1, Asia II 6, China 
1, MEAMl, and MED) of the 11 whitefly species from 
both laboratory and field. The data indicate that the 
frequencies of Wolbachia infection between laboratory 
and field populations in each of the five species appeared 



similar (Table Al). Thus, the laboratory rearing seemed to 
have exerted little effects on the frequencies of Wolbachia 
infection in these whitefly species. Until now, factors 
underlying the high variability of Wolbachia infection 
between the whitefly species are virtually unknown but 
certainly warrant future investigations. 

In contrast to a previous study that reports absence of 
Wolbachia infection in B. afer populations from China 
(Chu et al. 2010), the rate of Wolbachia infection in the 
B. afer population examined in the current study reached 
77.5%. Phylogenetic analysis of rrs, groEL, MLST, and wsp 
genes showed that the Wolbachia detected from B. afer 
belongs to supergroup B, which agrees with the report of 
Nirgianaki et al. (2003). 

Identification of a novel Wolbachia 
supergroup O 

Preliminary Bayesian phylogenetic analysis based on rrs 
gene sequences strongly supports the existence of one 
strange monophyletic group compared with the other 
Wolbachia identified in whiteflies. The rrs sequences from 
five of the whitefly populations {wBt_2, wBt_W, wBt_28- 
2, wBt_29-2, and wi?f_30-2) were clustered into group O. 
Average distance among those strange rrs sequences (592 bp) 
are 0.48%. The divergence of rrs between wBt_lO and all 
previously described Wolbachia supergroups (A to N) is 
higher than the 2% distance, a level of divergence that 
may merit the establishment of a new supergroup 
(Stouthamer et al. 1993; Augustinos et al. 2011). What is 
more, independent Bayesian analysis of rrs and groEL 
gene sequences and Neighbor-net analysis of gltA, hcpA, 
and fbpA gene sequences confirmed the distinct phyloge- 
netic position of M'5t_10 from the other supergroups. 
Based on the evidence, we propose the strange Wolbachia 
group as a new supergroup - Supergroup O. 
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(C) . I (D) 



Figure 7. Whole-mount FISH of Bemisia tabaci 
nymplns and female adults using a Portiera- 
specific probe (red) and a M/o/faach/a-specific 
probe (green). Upper column, Asia II 1 nymph 
and female adult; lower column, Asia II 9 
nymph and female adult. A, C, E, G: 
Wolbachia channel on a dark-field channel. B, 
D, F, H: Overlay of Portiera and Wolbachia 
channels on a bright-field channel. White 
triangles in D and H indicate anterior poles of 
the oocytes. Signals on legs, joints, and wings 
are chitin autofluoresence. 
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All previously known Wolbachia in Bemisia 
tabaci belong to supergroup B 

Except for the five supergroup O strains, phylogenetic 
analysis of eight molecular markers (rrs, groEL, gatB, 
coxA, hcpA, fbpA, ftsZ, and wsp genes) showed that all the 
Wolbachia strains detected from Chinese whiteflies as well 
as one strain from the Australia species belong to super- 
group B. This is consistent with previous diversity studies 
on Bemisia and Trialeurodes whiteflies (Nirgianaki et al. 
2003; Sintupachee et al. 2006; Gueguen et al. 2010; Singh 
et al. 2012; Tsagkarakou et al. 2012). 

The protein-coding genes are limited in 
Wolbachia diversity investigation 

At the early stage of Wolbachia research, the identification 
of Wolbachia strains was inferred based on the rrs gene 
(O'Neill et al. 1992; Stouthamer et al. 1993; Dumler and 
Walker 2005). As the research progressed, the rrs gene 
was found too conserved for further analysis of the Wol- 
bachia genus. Subsequently, additional protein-coding 
genes {gltA, groEL, ftsZ, and wsp genes) were developed 
for infection and evolutionary analysis of Wolbachia 
(Werren et al. 1995b; Zhou et al. 1998; Lo et al. 2002, 
2007; Casiraghi et al. 2005). Baldo et al. (2006) developed 
a standard MLST-based system (gatB, coxA, hcpA, ftsZ, 
and jhpA) for genotyping and strain classification of Wol- 
bachia infections. However, with more exploration of 
Wolbachia diversity, conflict results occurred among these 
different markers (Augustinos et al. 2011). In this study, 
the presence of supergroup O was confirmed by rrs and 
four protein-coding genes (fbpA, gltA, groEL, and hep 
genes) (Figs. 2-4, and Figs. AlO and A12). Whereas phy- 
logenetic analysis of several protein-coding genes {coxA, 
groEL, gatB, ftsZ, and wsp) clustered many Wolbachia O 
strains into supergroup B (Figs. 4 and 5, and Figs. A9, A8 
and All). Similar phenomena have been noticed in previ- 
ous studies. For example, even though the supergroup M 
and N have been identified as new groups of Wolbachia 
by rrs gene clustering, Augustinos et al. (2011) found that 
several popular protein-coding sequences such as gltA, 
groEL, and MLST genes clustered some individuals of 
those new groups into the old supergroup B. Besides, 
failures of amplifying MLST and wsp genes in many Wol- 
bachia O-infected whiteflies (Table 1) indicated protein- 
coding genes may not be sufficient for investigating the 
diversity of Wolbachia in B. tabaci. The failure of 
amplification of ftsZ and wsp genes were also observed in 
Wo/foac?i!fl-infected aphids (Augustinos et al. 2011). Con- 
sequently, it seems clear that phylogenetic analysis merely 
using protein-coding genes may underestimate the diver- 
sity of Wolbachia. 



That inadequacy of protein-coding genes for analyzing 
the diversity of Wolbachia may be explained by: (1) prim- 
ers of protein-coding genes are designed based on the ear- 
liest known Wolbachia (mostly A and B); and (2) different 
protein-coding genes suffer different selective pressure 
and thus have different evolutionary patterns. The rrs 
gene sequence is more conserved than wsp gene and the 
results of amplification are more stable compared with 
that of wsp or ftsZ genes which often produces unex- 
pected bands (not target-size bands or not the gene of 
Wolbachia). In fact, no single pair of primers can ensure 
detection of all Wolbachia specifically among various sam- 
ples (SimOEs et al. 2011). In view of the limitation of the 
various primers, we suggest that infection data obtained 
by any of these genes should be confirmed by vector clon- 
ing and sequencing of all representative bands. 

Wolbachia in Bemisia tabaci are transmitted 
horizontally 

Our FISH data indicate that Wolbachia can be vertically 
transmitted in whiteflies (Fig. 7), a result in agreement 
with that of a previous report (Gottlieb et al. 2008). In 
addition, our FISH data show the distribution of Wolba- 
chia outside of bacteriocytes of the whitefly adults and thus 
also indicate potential horizontal transmission of Wolba- 
chia. Not surprisingly, incongruence was found between 
the phylogeny of Bemisia mtCOI sequences and that of 
Wolbachia supergroup B based on concatenated MLST 
sequences (Fig. 6). In addition, in several cases, a popula- 
tion of a given whitefly species harbored divergent Wolba- 
chia strains (e.g., Population of 28 in Table Al). As 
speculated by a rate of rrs gene divergence of 1-2% per 50 
million years in bacterial endosymbionts (Moran et al. 
1993; Ochman et al. 1999), the divergence between super- 
group B and supergroup O probably started more than 
120 million years ago. While the divergence date of B. tab- 
aci complex was speculated to start about 50 million years 
ago, much more recent than that of supergroup B and O 
Wolbachia (Boykin et al. 2013). The double infection of 
Wolbachia supergroups B and O in the same population 
indicates horizontal transmission of Wolbachia. Horizontal 
transmission of Wolbachia has often been speculated based 
on phylogenetic analysis (Werren et al. 1995a; Sintupachee 
et al. 2006; Stahlhut et al. 2010; Schuler et al. 2013; Zhang 
et al. 2013). Wolbachia has also been reported from other 
whitefly genera such as Trialeurodes and some parasitoids 
(Raychoudhury et al. 2009; Cass et al. 2014), and this 
diversity of distribution may also hint horizontal transmis- 
sion. Sintupachee et al. (2006) hypothesized that the 
horizontal transmission of Wolbachia from whiteflies to 
other arthropods may occur through plants, because 
whiteflies could feed on plants without ruining plant cells. 
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Caspi-Fluger et al. (2011) presented a case study of hori- 
zontal transmission of Rickttesia in whiteflies via plants. 
Though we are yet unable to speculate on the origin of 
Wolbachia in whiteflies, we suggest that horizontal trans- 
mission of Wolbachia in whiteflies via plants warrants 
investigation especially as this bacterium has been detected 
outside of the bacteriocytes in the insect hosts. 

Conclusion 

We conducted a comprehensive screening for Wolbachia 
in whiteflies, and the findings have broadened substan- 
tially the host spectrum of Wolbachia and revealed a new 
supergroup of Wolbachia in whiteflies. Our study also 
shows the limitations of protein-coding genes as molecu- 
lar markers for Wolbachia investigation. Both specific and 
efficient molecular markers are needed for intensive sur- 
veys of Wolbachia. Wolbachia are transmitted vertically 
and horizontally in whiteflies. Clarifying the Wolbachia 
strains of whiteflies and their biological functions may 
provide novel clues for the development of efficient con- 
trol technologies against invasive whiteflies and whitefly- 
transmitted plant viruses. 
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Appendix 



Table A1. Details of screen of Bemisia afer and B. tabaci cryptic species for Wolbachia. Rates of infections do not differ between sexes in eacli 
of tine populations analyzed with Fisher's exact test (Populations with fewer than 10 individuals were excluded from the analysis). 





Whitefly 
species 








Collection 




Sample size 


Infection rate (%)' 




No. 


Location 


Latitude 


Longitude 


date 


Host plant (family)' 


F 


M UN 


Females 


Males 


Overall 


1 


Bemisia afer 


Linyi, Shandong, 
China 


35°47'N 


118°37'E 


July 2012 


6rau55onef/a 
papyfifera (1) 


24 


16 - 


87.50 


62.50 


77.50 


Bemisia 
























tabaci 
























2^ 


MED 


Hefei, Anhui, China 


3r95'N 


117°48'E 


October 2009 


Solanum 
melongena (2) 


21 


8 - 


9.52 


0.00 


6.90 


3 


MED 


Hefei, Anhui, China 


3r95'N 


117°48'E 


October 2009 


Solanum 
lycopersicum (2) 


1 7 
1 / 


1 J — 


n nn 
u.uu 


u.uu 


0.00 


4 


MEAM1 


Hefei, Anhui, China 


3r92'N 


117°14'E 


October 2009 


Salvia splendens (3) 


21 


7 - 


0.00 


0.00 


0.00 


5 


MED 


Nanjing, Jiangsu, 
China 


32°02'N 


118°B4'E 


September 
2009 


Solanum 
melongena (2) 


25 


5 - 


0.00 


0.00 


0.00 


6 


MED 


Nanjing, Jiangsu, 
China 


32°02'N 


118°54'E 


September 
2009 


Brassica oleracea 
var. capitata (4) 


15 


15 - 


0.00 


0.00 


0.00 


7 


MEAM1 


Chongmingdao, 
Shanghai, China 


3r50'N 


121°80'E 


November 
2009 


Brassica oleracea 
var. capitata (4) 


16 


11 - 


0.00 


0.00 


0.00 


8 


MED 


Chongmingdao, 
Shanghai, China 


3r50'N 


121°80'E 


November 
2009 


Capsicum annuum (2) 


21 


10 - 


0.00 


0.00 


0.00 


9" 


Asia II 3 


Hangzhou, Zhejiang, 
China 


30°23'N 


120°18'E 


April 2009 


Glycine max (5) 


26 


21 - 


0.00 


0.00 


0.00 


10^'" 


Asia II 1 


Hangzhou, Zhejiang, 
China 


29°27'N 


119°18'E 


October 2010 


Gossypium 
hirsutum (6) 


29 


15 - 


100.00 


100.00 


100.00 


ll'* 


China 1 


Hangzhou, Zhejiang, 
China 


30°23'N 


120°18'E 


November 
2009 


Solanum 
lycopersicum (2) 


23 


29 - 


100.00 


82.76 


90.38 


12'' 


MED 


Ningbo, Zhejiang, 
China 


29°48'N 


121°3B'E 


June 2009 


Capsicum annuum (2) 


22 


16 - 


0.00 


0.00 


0.00 


13 


MED 


Taizhou, Zhejiang, 
China 


28°30'N 


121°34'E 


October 2012 


Cucurbita moschata (7) 


22 


11 - 


31.82 


9.09 


24.24 


14^ 


MEAM1 


Wenzhou, Zhejiang, 
China 


27°47'N 


120°39'E 


September 
2008 


Solanum melongena (2) 


20 


20 - 


0.00 


0.00 


0.00 


15 


MED 


Nanchang, Jiangxi, 
China 


28°72'N 


115°91'E 


October 2009 


Cucurbita moschata (7) 


16 


11 - 


0.00 


0.00 


0.00 


16 


MED 


Nanchang, Jiangxi, 
China 


28°72'N 


115°91'E 


October 2009 


Ipomoea batatas (8) 


16 


14 - 


0.00 


0.00 


0.00 


17 


MED 


Nanchang, Jiangxi, 
China 


28°72'N 


115°91'E 


October 2009 


Brassica campestris 
ssp. Pekinensis (4) 


21 


8 - 


0.00 


0.00 


0.00 



2728 



© 2014 The Authors. Ecology and Evolution published by John Wiley Si Sons Ltd. 



X. -L. Bing ef al. Wolbachia in Whiteflies 



Table Al. Continued. 



No. 


Whitefly 
species 


Location 


Latitude 


Longitude 


Collection 
date 


Host plant (family)' 


Sample size 
F M UN 


Infection rate (%)■ 
Females Males 


Overall 


18 


MED 


Nanchang, Jiangxi, 


28°72'N 


115°91'E 


October 2009 


Humulus scandens (9) 


16 


12 - 


0.00 


0.00 


0.00 






China 




















19 


MED 


Nanchang, Jiangxi, 


28°29'N 


116°01'E 


October 2009 


Citrullus lanatus (7) 


13 


16 - 


0.00 


0.00 


0.00 






China 




















20 


MED 


Nanchang, Jiangxi, 


28°29'N 


116°01'E 


October 2009 


Ipomoea batatas (8) 


20 


7 - 


0.00 


0.00 


0.00 






China 




















21 


MED 


Jiujiang, Jiangxi, 


29°76'N 


115°79'E 


October 2009 


Capsicum annuum (2) 


16 


11 - 


0.00 


0.00 


0.00 






China 




















22 


MED 


Jiujiang, Jiangxi, 


29°76'N 


115°79'E 


October 2009 


Ipomoea batatas (8) 


24 


6 - 


0.00 


0.00 


0.00 






China 




















23 


MED 


Jiujiang, Jiangxi, 


29°76'N 


115°79'E 


October 2009 


Cucumis sativus (7) 


20 


8 - 


0.00 


0.00 


0.00 






China 




















24 


MED 


Jiujiang, Jiangxi, 


29°76'N 


115°79'E 


October 2009 


Solanum melongena (2) 


13 


2 - 


0.00 


0.00 


0.00 






China 




















25 


MED 


Jiujiang, Jiangxi, 


29°76'N 


115°79'E 


October 2009 


Phaseolus vulgaris (5) 


16 


11 - 


0.00 


0.00 


0.00 






China 




















26" 


Asia II 7 


Guangzhou, 


23°09'N 


113°21'E 


October 2007 


Gossypium hirsutum (5) 


30 


8 - 


100.00 


87.50 


97.37 






Guangdong, 
























China 




















27" 


Asia 1 


Zhaoqing, 


23°56'N 


112°1'E 


November 


Ipomoea batatas (8) 


20 


20 - 


100.00 


100.00 


100.00 






Guangdong, 






2010 


















China 




















28^ 


Asia II 1 


Zhaoqing, 


23°56'N 


112°1'E 


August 2012 


Arachis hypogaea (5) 


31 


12 - 


96.77 


83.33 


93.02 






Guangdong, 
























China 




















29^ 


Asia II 1 


Zhaoqing, 


23°56'N 


112°1'E 


August 2012 


Ipomoea batatas (8) 


32 


3 - 


96.88 


100.00 


97.14 






Guangdong, 
























China 




















30^ 


Asia II 1 


Sanya, Hainan, 


18°24'N 


109°42'E 


August 2009 


Ipomoea batatas (8) 


24 


5 - 


100.00 


80.00 


96.55 






China 




















31 


MEAM1 


Nanning, Guangxi, 


22°38'N 


108°23'E 


August 2009 


Vigna unguiculata (5) 


17 


13 - 


0.00 


0.00 


0.00 






China 




















32 


MEAM1 


Nanning, Guangxi, 


22°38'N 


108°23'E 


August 2009 


Gossypium 


22 


0 - 


0.00 




0.00 






China 








hirsutum (6) 












33 


MEAM1 


Nanning, Guangxi, 


22°38'N 


108°23'E 


August 2009 


Cucumis sativus (7) 


17 


13 - 


0.00 


0.00 


0.00 






China 




















34 


MED 


Beihai, Guangxi, 


21°29'N 


109°09'E 


August 2009 


Ipomoea batatas (8) 


10 


12 - 


0.00 


0.00 


0.00 






China 




















35 


Asia II 6 


Baise, Guangxi, 


22°94'N 


108°54'E 


August 2009 


Luffa cylindrica (7) 


15 


3 - 


53.33 


0.00 


44.44 






China 




















36 


MEAM1 


Baise, Guangxi, 


23°52'N 


106°37'E 


August 2009 


Vigna unguiculata (5) 


4 


9 - 


0.00 


0.00 


0.00 






China 




















37 


MEAM1 


Baise, Guangxi, 


23°45'N 


106°47'E 


August 2009 


Benincasa hispida (7) 


15 


12 - 


0.00 


0.00 


0.00 






China 




















38" 


Asia II 6 


Baise, Guangxi, 


22°94'N 


108°54'E 


November 


Ipomoea batatas (8) 


36 


6 - 


13.89 


16.67 


14.29 






China 






2011 
















Asia 1 


Honghe, Yunnan, 


z4 JO IM 


1 Ud 4d b 


November 


Ipomoea batatas (8) 


46 


37 - 


97.83 


100.00 


QQ QD 

yo.oU 






China 






2011 














40 


MED 


Guiyang, Guizhou, 


26°40'N 


106°67'E 


July 2009 


Glechoma longituba (3) 


9 


21 - 


0.00 


0.00 


0.00 






China 




















41 


China 1 


Zunyi, Guizhou, 


27°39'N 


107°70'E 


July 2009 


Solanum melongena (2) 


24 


6 - 


75.00 


83.33 


76.67 






China 




















42 


China 1 


Zunyi, Guizhou, 


27°39'N 


107°70'E 


July 2009 


Ipomoea batatas (8) 


23 


7 - 


100.00 


85.71 


96.67 






China 
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No. 


Winitefly 
species 


Location 


Latitude 


Longitude 


Collection 
date 


Host plant (family)' 


F 


M UN 


Females 


Males 


Overall 


43^ 


Asia II 9 


Shaoyang, Hunan, 


26°59'N 


111°16'E 


October 201 1 


Ipomoea batatas (8) 


19 


16 - 


89.47 


87.50 


88.57 






Cliina 




















44 


MED 


Jishou, Hunan, 


28°18'N 


109°38'E 


September 


Raphanus sativus (4) 


7 


8 - 


0.00 


0.00 


0.00 






China 






2009 














45 


MED 


Luoyang, Henan, 


34°46'N 


112°46'E 


September 


Ipomoea batatas (8) 


16 


5 - 


0.00 


0.00 


0.00 






Cliina 






2009 














46 


MEAM1 


Luoyang, Henan, 


34°66'N 


112°51'E 


September 


Gossypium 


15 


15 - 


0.00 


0.00 


0.00 






China 






2009 


hirsutum (6) 












47 


MEAM1 


Luoyang, Henan, 


34°59'N 


112°58'E 


September 


Solan urn 


20 


8 - 


0.00 


0.00 


0.00 






China 






2009 


melongena (2) 












48 


MEAM1 


Luoyang, Henan, 


34°59'N 


112°58'E 


September 


Cucumis sativus (7) 


26 


4 - 


0.00 


0.00 


0.00 






China 






2009 














49-1 


Asia II 3 


Zhengzhou, Henan, 


34°78'N 


113°66'E 


September 


Solanum melongena (2) 


0 


1 - 


- 


100.00 


100.00 






China 






2009 














49-2 


MED 


Zhengzhou, Henan, 


34°78'N 


113°66'E 


September 


Solanum melongena (2) 


18 


9 - 


0.00 


0.00 


0.00 






China 






2009 














50 


MEAM1 


Zhengzhou, Henan, 


34°78'N 


113°66'E 


September 


Cucumis sativus (7) 


27 


2 - 


0.00 


0.00 


0.00 






China 






2009 














51 


MED 


Xinxiang, Henan, 


35°47'N 


1 13°75'E 


September 


Phaseolus vulgaris (5) 


12 


4 - 


0.00 


0.00 


0.00 






China 






2009 














52 


MED 


Xinxiang, Henan, 


35°47'N 


113°75'E 


September 


Raphanus sativus (4) 


26 


3 - 


0.00 


0.00 


0.00 






China 






2009 














53 


MED 


Xinxiang, Henan, 


35°47'N 


113°75'E 


September 


Brassica campestris 


20 


8 - 


0.00 


0.00 


0.00 






China 






2009 


ssp. Pekinensis (4) 












54 


Australia 


Bundaberg, 


24°48'S 


152°27'E 




Euphorbia 




- 11 






100.00 






Queensland, 








cyathophora (10) 
















Australia 





















F, female adult; M, male adult; UN, unknown sex; -, not ascertained. 

'in all 22 species of host plants from 10 families, figures in parentheses indicate the names of the families: (1), Moraceae; (2), Solanaceae, (3), 
Lamiaceae, (4), Cruciferae, (5), Fabaceae, (6), Malvaceae, (7), Cucurbitaceae, (8), Convolvulaceae, (9), Cannabaceae, (10), Euphorbiaceae. 
^Infection rates of Wolbachia detected by diagnostic PCR of rrs gene. 
^Populations for the detection of Wolbachia supergroup 0 (Table 3). 
''Populations maintained in the laboratory on cotton since collection. 



Table A2. List of the primers used for screening and sequencing. 



Gene 



Hypothetical product 



Primer name 



Primer sequences (5'-3') 



Tm Product size Reference 



Bemlsia spp. 

mtCOl Mitochondrial cytochrome 
oxidase subunit I 
Universal bacteria 

rrs Ribosomal RNA 16S 

Wolbachia spp. 
rrs Ribosomal RNA 16S 

gatB Glutamyl-tRNA(Gln) 

amidotransferase, subunit B 
coxA Cytochrome coxidase, subunit I 

hcpA Conserved hypothetical protein 



COI-F-C1-J-2195: TTGATninGGTCATCCAGAAGT 54°C 759 bp 
COI-R-TL2-N-3014: TCCAATGCACTAATCTGCCATATTA 



27F: 
1494R: 

W0I-I6S-F: 

W0I-I6S-R: 

gatB_F1: 

gatB_R1: 

coxA_F1 : 

coxA_R1: 

hcpA_F1: 

hcpA_R1: 



AG AGTHG ATC MTG G CTC AG 
CTAC G G CTAC CTTGnAC G A 

CGGGGGAAAAATTTATTGCT 

AG CTGTAATAC AG AAAGTAAA 

GAKTTAAAYCGYGCAGGBGn 

TGGYAAYTCRGGYAAAGATGA 

TTGGRGCRATYAACTTTATAG 

CTAAAGACHTKACRCCAGT 

GAAATARCAGTTGCTGCAAA 

GAAAGTYRAGCAAGYTCTG 



55°C 589 bp 

54°C 471 bp 

54°C 487 bp 

54°C 515 bp 



Frohlich et al. (1999) 



50°C 1417 bp Weisburg et al. (1991) 



Heddi et al. (1999) 
Baldo et al. (2006) 
Baldo et al. (2006) 
Baldo et al. (2006) 
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Gene 


Hypothetical product 


Primer name 


Primer sequences (B'-3') 


Tm 


Product size 


Reference 


ftsZ 


Cell division protein 


ftsZ_F1: 


ATYATGGARCATATAAARGATAG 


54°C 


524 bp 


Baldo et al. (2006) 






ftsZ_R1: 


TC RAG YAATG G ATTRG ATAT 








fbpA 


rr ULLUitr uij|Ji ius[Jr laitr aiuuiajtr 


fbpA Fl ■ 




59°C 


509 bp 








fbpA_R1: 


CCRCCAGARAAAAYYACTATTC 








wsp 


Outer surface protein 


wsp_F1 : 


GTCCAATARSTGATGARGAAAC 


B9°C 


546 bp 


Baldo et al. (2006) 






wsp_R1 : 


CYGCACCAAYAGYRCTRTAAA 








groEL 


Chaperonin GroEL 


groEL-F: 


CAACRGTRGSRRYAACTGCDGG 


B4°C 


491 bp 


Ros et al. (2009) 






groEL-R: 


GATADCCRCGRTCAAAYTGC 








gItA 


Citrate synthase 


WgltAFI: 


TACGATCCAGGGTnGTTTCTAC 


B4°C 


659 bp 


Casiraghi et al. (2005) 






WgltARev2: 


CATTTCATACCACTGGGC 









0.96 



1.00 



0.65 



0.98 



0.98 



0.97 



1.00 



Pop. 30-2 (KF452521) 
Pop. 30-1 {KF452522) 
Pop. 28-1 (KF452517) 
Asia II 1 (DQ309077) 
Pop. 10 (KF452529) 
Pop. 29-1 (KF452520) 
Pop. 29-2 (KF452519) 
Pop. 28-2 (KF452518) 
Q g2p Asia II 6 (AJ784261) 
Pop. 38 (KC540758) 



1.00 



0.89 



0.83 



Pop. 35 (KF452523) 

Asia II 5 (DQ1 16673) 
Pop. 26 (EUig2043) 
1-001 Asia II 7 {AJ748378) 
Asia II 2 (AY686088) 
QQ, Pop. 49-1 (KF452527) 
T- Asia II 3 (AJ867556) 
1:00 Pop. 43 (KF452530) 
I.OOI Asia II 9 (HM137313) 
Pop. 41 (KF452524) 
China 1 (AY686089) 
Pop. 11 (GQ303180) 
Pop. 42 (KF452525) 
China 2 (AY686072) 
Pop. 27 (KF452531) 
Pop. 39 {KC540757) 
Asia I (AJ748370) 



1.00 
1.00 



1.00 



0.99 



jjjQT Pop. 54 (KF452532) 

' I 



1.00 
1.00 



- Australia (DQ1 30052) 
Pop. 13 (KF452516) 
■ Mediterranean (FJ766381) 

Pop. 2 (KF452526) 
IVliddle East-Asia IVIInor 1 (AF418671) 
Pop. 1 (KF452528) 
^ Bemisia afer (GU220055) 
Trialeurodes vaporariorum (AY521265) 



0.1 



Figure A3. Phylogenetic analysis of the Bemisia spp. based on whitefly mtCOl gene sequences (657 sites). Trialeurodes vaporariorum is used as out 
group. Reference sequences are obtained from the CSIRO data access portal (De Barro and Boykin 2013). The tree was constructed using a 
TIM3 + I + G substitution model for Bayesian analysis. Bayesian posterior probabilities are indicated at nodes. The sequences obtained in this study 
are shown in bold. The bar indicates a branch length of 0.1 substitutions/site. The sequence GenBank accession numbers are shown in parentheses. 
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Phylum 


Class 


Order 


Host species 


15S rRNA gene 


gItA 


groEL 


Supergroup 


Arthropoda 


Insecta 


Hymenoptera 


Muscidifurax uniraptor 


L02882 






A 


Arthropoda 


Insecta 


Hymenoptera 


Nasonia vitripennis 


M84688 


AY7 14795 


AY714812 


A 


Arthropoda 


Prostigmata 


Acarina 


Bryobia 5arothamni 


EU499315 




EU499330 


B 


Arthropoda 


Prostigmata 


Acarina 


Bryobia praetiosa 


EU499317 


EU499327 


EU499332 


B 


Arthropoda 


Insecta 


Hymenoptera 


Nasonia vitripennis 


M84686 


AY7 14782 


AY7 14796 


B 


Arthropoda 


Insecta 


Hemiptera 


Bemisia tabaci 


JN204507 






B 


Nematoda 


Secernentea 


Spirurida 


Onchocerca ochengi 


AJ0 10276 


AJ609640 




C 


Nematoda 


Secernentea 


Spirurida 


Onchocerca gibsoni 


AJ276499 


AJ609639 


AJ6096B2 


C 


Nematoda 


Secernentea 


Spirurida 


Dirofilaria repens 


AJ276B00 




AJ609653 


c 


Nematoda 


Secernentea 


Spirurida 


Dirofilaria immitis 


Z49261 


AJ609641 




c 


Nematoda 


Chromadorea 


Spirurida 


Brugia malayi 


AF05114B 


AJ 609643 


AE017321 


D 


Nematoda 


Secernentea 


Spirurida 


Litomosoides sigmodontis 


AF069068 


AJ609645 


AF409113 


D 


Arthropoda 


Collembola 


Collembola 


Foisomia Candida 


API 79630 


AJ609649 




E 


Arthropoda 


Ellipura 


Collembola 


Mesaphorura macrocheta 


AJ422184 






E 


Arthropoda 


Insecta 


Neuroptera 


Myrmeleon mobiiis 


DQ068882 






F 


Arthropoda 


Insecta 


Isoptera 


Kalotermes flavicollis 


Y11377 


AJ609651 


AJ609660 


F 


Nematoda 


Secernentea 


Spirurida 


Ivlansonelia ozzardi 


AJ279034 




AJ6096B7 


F 


Arthropoda 


Insecta 


Isoptera 


Zootermopsis nevadensis 


AY764280 


AY764282 


AY764277 


H 


Arthropoda 


Insecta 


Siphonaptera 


Ctenocephaiides feiis 


AY335923 


AJ609650 


AJ6096B9 


1 


Arthropoda 


Insecta 


Siphonaptera 


Orchopeas ieucopus 


AY335924 






1 


Nematoda 


Secernentea 


Spirurida 


Dipetalonema graciie 


AJ548802 


AJ609648 


AJ6096B8 


J 


Arthropoda 


Arachnida 


Prostigmata 


Bryobia sp. 


EU499316 


EU499326 


EU499331 


K 


Nematoda 


Phasmida 


Tylenchida 


Radophoius similis 


EU833482 




EU833484 


L 


Arthropoda 


Insecta 


Hemiptera 


Tuberoiachnus saiignu 


JN384085 






M 


Arthropoda 


Insecta 


Hemiptera 


Aphis sp. 


JN384091 






M 


Arthropoda 


Insecta 


Hemiptera 


Cinara cedri 






JN3840B3 


M 


Arthropoda 


Insecta 


Hemiptera 


Toxoptera aurantii 


JN384094 






N 


Arthropoda 


Insecta 


Hemiptera 


Toxoptera aurantii 


JN384095 






N 


Arthropoda 


Insecta 


Hemiptera 


Bemisia tabaci 


KF454771 


KF587270 


KF452543 


O 




Figure A5. Phylogenetic position of the 
Wolbachia identified from Bemisia tabaci 
putative species Asia II 1 {wBt_^ 0) based on 
bacterial rrs gene sequences (1317 sites). 
Wolbachia strains are characterized by the 
names of their host species. The tree was 
constructed using a HKY + G substitution 
model for Bayesian analysis. Bayesian posterior 
probabilities are shown on the branches. The 
sequence obtained in this study is shown in 
bold. The bar indicates a branch length of 0.01 
substitutions/site. The names and sequence 
GenBank accession numbers are listed in 



Tables A4 and A6. 
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M 1 2 3 4 5 6 M 8 9 10 11 12 13 



Table A6. Divergence between the Wolbachia detected in the puta- 
tive species Asia II 1 of the 6. tabaci complex (wBtJlO) and other 
supergroups based on 16S rRNA gene sequences. 

GenBank 



Host species 


Supergroup 


Divergence {%) 


accession no. 


Muscidifurax uniraptor 


A 


3.73 


L02882 


Nasonia vitripennis 


A 


3.64 


M84688 


Bryobia sarothamni 


B 


4.98 


EU499315 


Bryobia praetiosa 


B 


4.98 


EU499317 


Bemisia tabaci 


B 


5.18 


JN204507 


Nasonia vitfipennis 


B 


4.83 


M84686 


Onchocerca ochengi 


C 


5.61 


AJ010276 


Onchocerca gibsoni 


C 


5.44 


AJ276499 


Dirofilaria repens 


C 


5.99 


AJ276500 


Dirofilaria immitis 


C 


5.63 


Z49261 


Brugia malayi 


D 


4.59 


AF051145 


Litomosoides sigmodontis 


D 


5.15 


AF069068 


Foisomia Candida 


E 


3.65 


AF1 79630 


l\/Iesapi)orura macrocheta 


E 


4.04 


AJ422184 


IVIansonelia ozzardi 


F 


4.95 


AJ279034 


Myrnneleon mobilis 


F 


3.90 


DQ068882 


Kaiotermes flavicollis 


F 


3.70 


Y11377 


Zootermopsis nevadensis 


H 


4.41 


AY7 64280 


Ctenocephalides feiis 


1 


7.22 


AY335923 


Orchopeas ieucopus 


1 


7.48 


AY335924 


Dipetalonema graciie 


J 


5.59 


AJ548802 


Bryobia sp. 


K 


3.25 


EU499316 


Radopholus similis. 


L 


4.79 


EU833482 


Tuberolachnus salignu 


M 


2.52 


JN384085 


Aphis sp. 


M 


2.59 


JN384091 


Toxoptera aurantii 


N 


4.50 


JN384094 


Toxoptera aurantii 


N 


4.80 


JN384095 



Figure A7. RFLP pattern of PGR products of rrs gene of the 
Woibachia supergroup 0 and B in Bemisia tabaci corresponding to 
Vspl (Asel) digestion. The different profiles were obtained from 
individuals representing different Wolbachia in 6. tabaci. The bands 
shown on the bottom are primer dimers. Lane 1, undigested 
Wolbachia 0 obtained by PGR from Asia II 1; lane 2-4, undigested 
Wolbachia B obtained by PGR from Asia II 7, Asia I and Ghina 1, 
respectively; lane 5-6, Wolbachia PGR amplified production from 
MEAM1 and IVIED as controls; lane 8, digested Wolbachia 0 obtained 
by PGR from Asia II 1; lane 9-11, Wolbachia B obtained by PGR from 
Asia II 7, Asia I and Ghina 1, respectively; lane 12-13, Wolbachia PGR 
amplified production from MEAMl and MED as controls; lane M, 
DNA size markers (100, 250, 500, 750, 1000, and 2000 bp from 
bottom to top). 
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1.00 



Brugia malayi (28) 

Zootermopsis angusticollis (64) 

Cimex lectularius (26) 



D 
H 



0.95 



0.1 



1.00 



- Opistophthalmus granifrons (30) 

i- Glossina morsitans gambiensis (145) 
- Bemisia tabaci (9) 
-wBf_27 (210) 
-wBaJ, wBt_^^, 26 (105) 
-wBf_35 (216) 
- wSt_30-1 (208) 

-wSMO, 28-1, 29-1, 38, 39, 41, 43, 49-1, 54 (207) 
-wBf_42 (209) 
-Pheidole micula (72) 



0.73!^ 



0.75 



■'V 

0.95L 



1.00 



— Glossina morsitans morsitans (141) 

Drosoptiila simulans (22) 

L Drosophila melanogaster (1 ) 

Figure A8. Phylogenetic position of tine Wolbactiia identified from Bemisia afer and 6. tabaci putative species based on bacterial gatB gene 
sequences (369 sites). Woibacliia strains are cinaracterized by tine names of their liost species and allele numbers from the MLST database. The 
tree was constructed using a TPM2uf + G substitution model for Bayesian analysis. Bayesian posterior probabilities are shown on the branches. 
The sequences obtained in this study are shown in bold. The bar indicates a branch length of 0.1 substitutions/site. 



Brugia malayi (29) 

Zootermopsis angusticollis (54) 
Drosophila melanogaster (1 ) 
Glossina morsitans morsitans (23) 



D 
H 



1.00 



Opi 



1.00 



°-^'^LCaiyptrate sp. (47) 

Cimex lectularius (27) 
Opistoptithaimus granifrons (31 ) 
iv8a_1 (11) 

_ Teleogrylius talwanemma (25) 
Tribolium confusum (5) 

26, 27, 28-1, 28-2, 29-1, 30-1, 30-2, 35, 38, 39, 41, 42, 43, 49-1, 54 ( 



0.1 



_ Tribolium 

1.00 

LwBt 11, 



Figure A9. Phylogenetic position of the Woibachia identified from Bemisia afer and B. tabaci putative species based on bacterial coxA gene 
sequences (402 sites). Woibachia strains are characterized by the names of their host species and allele numbers from the MLST database. The 
tree was constructed using a TIM1 + G substitution model for Bayesian analysis. Bayesian posterior probabilities are shown on the branches. The 
sequences obtained in this study are shown in bold. The bar indicates a branch length of 0.1 substitutions/site. 
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Figure A10. Phylogenetic position of the Wolbachia identified from Bemisia afer and 6. tabaci putative species based on bacterial hcpA gene 
sequences (444 sites) using tine Neighbor-net method. Each edge (or a set of parallel edges) corresponds to a split in the data set and has length 
equal to the weight of the split. The sequence obtained in this study is shown in bold. MLST Database allele numbers of hcpA sequences are 
shown in parenthesis. The bar indicates a branch length of 0.01 substitutions/site. 



0.87 



Brugia malayi (26) 

- Cimex lectularius (24) 

- Opistophthalmus granifrons (18) 
Zootermopsis angusticollis (57) 



0.82 



1.00 



j-Pheidole micula (17) 



0.95 



0.85 



0.1 



Drosophila melanogaster (1 ) 
.Drosophila simulans (3) 
0.90L Glossina morsitans morsitans (114) 

wBMO, 11, 28-1, 29-1, 30-1, 38, 41,42, 43, 54 (170) 
_wSf_49-1 (180) 
_wBf_39 (182) 
- Tribolium confusum (18) 
.wBM3 (181) 
wet_26, 27, 35 (7) 
Teleogryllus taiwanemma (20) 
Glossina morsitans gambiensis (117) 



1.00 



0.67 



O.gstwBaJ (11) 

Figure All. Phylogenetic position of the Wolbachia identified from Bemisia afer and B. tabaci putative species based on bacterial ftsZ gene 
sequences (435 sites). Wolbachia strains are characterized by the names of their host species and allele numbers from the MLST database. The 
tree was constructed using a TrN + I substitution model for Bayesian analysis. Bayesian posterior probabilities are shown on the branches. The 
sequences obtained in this study are shown in bold. The bar indicates a branch length of 0.1 substitutions/site. 
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Figure A12. Phylogenetic position of tine Wolbachia identified from Bemisia afer and 6. tabaci putative species based on bacteriai fbpA gene 
sequences (429 sites) using tlie Neiglnbor-net method. Eacii edge (or a set of paraiiei edges) corresponds to a spiit in tine data set and Inas iength 
equai to tine weigint of the spiit. The sequences obtained in this study are shown in boid. MLST Database aiieie numbers of fbpA sequences are 
shown in parenthesis. The bar indicates a branch iength of 0.01 substitutions/site. 
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Asia II 1 with no probe Asia II 1 without Wolbachia probe Asia II 1 without Wolbachia probe 



(G) (H) I (I) 




Mediterranean Milddle East-Asia Minor 1 Milddle East-Asia Minor 1 



Figure A13. FISH controls. No probe control: Asia II 1; No Wolbachia probe control: Asia II 1; Wolbachia-iree whiteflies control: Mediterranean 
(harboring Portiera and Hamiltonella) and Middle East-Asia Minor 1 (harboring Portiera, Hamiltonella and Rickettsia; Bing et al. 2013a). A-C, G-l: 
Overlay of channels of Portiera (red) and Wolbacliia (green); D-F, J-L: Overlay of channels of Portiera (red), Wolbachia (green) and white light. 
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